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Abstract 

The S 2 x R geometry can be derived by the direct product of the spher- 
ical plane S 2 and the real line R. In (H J. Z. Farkas has classified and 
given the complete list of the space groups of S 2 xR. The S 2 xR manifolds 
were classified by E. MOLNAR and J. Z. FARKAS in ||2l by similarity and 
diffeomorphism. In Ifl2l we have studied the geodesic balls and their vol- 
umes in S 2 x R space, moreover we have introduced the notion of geodesic 
ball packing and its density and have determined the densest geodesic ball 
packing for generalized Coxeter space groups of S 2 x R. 

In this paper we study the locally optimal ball packings to the S 2 x R 
space groups having Coxeter point groups and at least one of the generators 
is a glide reflection. We determine the densest simply transitive geodesic 
ball arrangements for the above space groups, moreover we compute their 
optimal densities and radii. 

The density of the densest packing is m 0.80407553, may be surpris- 
ing enough in comparison with the Euclidean result -5= ss 0.74048. E. 



*AMS Classification 2010: 52C17, 52C22, 53A35, 51M20 
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MOLNAR has shown in |4l, that the homogeneous 3-spaces have a unified 
interpretation in the real projective 3-sphere P5 3 (V 4 , V^M). In our work 
we shall use this projective model of S 2 x R geometry. 

1 Introduction 

S 2 x R is derived as the direct product of the spherical plane S 2 and the real line 
R. The points are described by (P, p) where P E S 2 and p E R. The isometry 
group 7som(S 2 xR) of S 2 xR can be derived by the direct product of the isometry 
group of the sphere Isom(S 2 ) and the isometry group of the real line Jsom(R). 

Isom(S 2 ) := {A E 0(3) : S 2 ^ S 2 : (P,p) i— >■ (PA,p)} for any fixed p. 

Isom(R) := {p : (P,p) (P, ±p + r)}, for any fixed P. 

r 

here the "-" sign provides a reflection in the point - E R, 
by the "+" sign we get a translation of R. 

(1.1) 

The structure of an isometry group T C Jsom(S 2 x R) is the following: T := 
{(Ax x pi),...(A n x p n )}, where Aj x ^ := At x (P^n) := (gi,n), (i E 
{1, 2, . . . n} and E Lsom(S 2 ), P« is either the identity map 1 R of R or the 
point reflection 1 R . g t := A^ x Pj is called the linear part of the transformation 
(Ai x and r, is its translation part. The multiplication formula is the following: 

(A, x P 1; n) o (A 2 x P 2 ,r 2 ) = ((A^a x R.R^r,^ + r 2 ). (1.2) 



Definition 1.1 A group of isometries Y C Jsom(S 2 x R) z'^ called space group 
if the linear parts form a finite group r called the point group of T, moreover, 
the translation parts to the identity of this point group are required to form a one 
dimensional lattice L r ofH. 

Remark 1.2 1. It can be proved that the space group T has a compact funda- 
mental domain T-p. 

2. IfT is not assumed to have a lattice, then it may have an infinite point group 

r . 
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Definition 1.3 The S 2 xR space groups T\ and T 2 are geometrically equivalent, 
called equivariant, if there is a "similarity" transformation X := S x o (S G 
Isom(S 2 ), a G S'zm(R)), such that T 2 = £ _1 ri£. Here a(s, t) : p — )■ p ■ s + t 
is a similarity of R, i. e. multiplication by 7^ s G R an J ?/zen addition by t G R 
/or every pGR 

Remark 1.4 T/' anJ T 2 are equivariant space groups then the their factor 
groups T\/Lr 1 and T 2 / L^ 2 are also equivariant. 

Thus the structure of the space group remains invariant under a similarity in the 
R-component and the spherical part is uniquely determined up to an isometry of 
S 2 . 

We characterize the spherical plane groups by the Macbeath- signature (see 

01, ma). 

In this paper we deal with such a S 2 x R space group where the generators 
gi, (i = 1,2, .. .m) of its point group T are reflections and at least one of the 
possible translation parts of the above generators unequal to zero. These groups 
are called glide reflection groups. 

Remark 1.5 In / f72l/ we have introduced the notion of generalized Coxeter group, 
if the generators gj, (i = 1,2, . . . m) of its point group T are reflections with 
translation parts t { = 0, (i = 1, 2, . . . m). 

In this paper we deal with the glide reflection space groups in S 2 x R space 
which are by denotation of [[I] : 

1. (+, 0, [ ]{(y)})xl a) q>2, 

r = (gi,g2 - gi,g|, (siS2) q ), 

2q.I.2: (f, I); 2qe. 1. 3: (0,1); 

2. (+, 0, [ ]{(2,2,g)})xl R , q>2, 

T = (gl, g2, g3 - g?, g|, gi, (gigs) 2 , (g2g 3 ) 2 ), (glg2) 9 ), 

4q.I.2: (0, 0, |); 4q. I. 3: (1, \, 0); 4q. 1. 4: {\,\,\); 
4qe.I.5: (0, \, 0); 4qe. 1. 6: (0,1,1); 

3. (+, 0, [] {(2,3,3)}) xl R , 

r = (gi, g 2 , g 3 - g 2 , gi, g 2 , (glg2) 2 , (gig 3 ) 3 , (g2g 3 ) 3 ), 
11. 1. 2: (|, |, |); 
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4. (+, 0, [ ] {(2,3,4)}) xl R , 

r = (gi, g 2 , g3 - gf, gi, gi, (gig2) 2 , (gig 3 ) 3 , (g2g 3 ) 4 ), 

12.1.2: (0, |,0); 12.1.3: (~, 0, ~); 12. 1. 4: (±,§,f); 

5. (+, 0, [ ] {(2,3,5)}) xl R , 

r = (gi, g2, g 3 - gf, gi, gf, (gig2) 2 , (gig 3 ) 3 , (g2g 3 ) 5 ), 

X. A. 2' 2/' 



2 Geodesic curve and balls in S 2 x R space 

E. Molnar has shown in 01, that the homogeneous 3-spaces have a unified inter- 
pretation in the projective 3-sphere VS 3 (V 4 , V4, R). In our work we shall use 
this projective model of S 2 x R and the Cartesian homogeneous coordinate sim- 
plex J B (e ), J E; i oo (e 1 ), J E; 2 3O (e2), £|°(e 3 ), ({e,} C V 4 with the unit point E(e = 
e o + e i + e 2 + e 3)) which is distinguished by an origin E and by the ideal 
points of coordinate axes, respectively. Moreover, y = cx with < c G I (or 
c G I \ {0}) defines a point (x) = (y) of the projective 3-sphere VS 3 (or that 
of the projective space V 3 where opposite rays (x) and (— x) are identified). The 
dual system {(e 4 )} C V 4 describes the simplex planes, especially the plane at 
infinity (e°) = E^E^Ef, and generally, v = u- defines a plane (u) = (v) of 
VS 3 (or that of V 3 ). Thus = xm = y v defines the incidence of point (x) = (y) 
and plane (it) = (v), as (x)I(w) also denotes it. Thus S 2 x R can be visualized in 
the affine 3-space A 3 (so in E 3 ) as well. 

In this context E. Molnar 01 has derived the well-known infinitezimal arc- 
length square at any point of S 2 x R as follows 

x 1 + y 2 + z 1 

We shall apply the usual geographical coordiantes (0, 6), (— it < 4> < tt, — | < 
9 < f ) of the sphere with the fibre coordinate i 6 R. We describe points in the 
above coordinate system in our model by the following equations: 

a; = 1, x 1 = e* cos cos 6 1 , x 2 = e* sin cos ^, x 3 = e t sin^. (2.2) 

Then we have x=^ = x 1 ,y = ^ = x 2 ,z = ^ J = x 3 , i.e. the usual Cartesian 
coordinates. We obtain by [4] that in this parametrization the infinitezimal arc- 
length square at any point of S 2 x R is the following 

{ds) 2 = (dt) 2 + (d<pf cos 2 e + (def. (2.3) 
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The geodesic curves of S 2 xR are generally defined as having locally minimal arc 
length between their any two (near enough) points. The equation systems of the 
parametrized geodesic curves j(t(r), (f>(r), 9(r)) in our model can be determined 
by the general theory of Riemann geometry (see [fi2l ). 

Then by (2.2) we get with c = sint>, u = cosv the equation systems of a 
geodesic curve, visualized in Fig. 1 in our Euclidean model: 



x(t) = e TSmv cos (r cosv), 

u, 

(2.4) 



y(r) = e TSmv sin (r cosv) cos u 



z(t) = e TSmv sin (r cosf ) smu 

71 71 
— 71 < U < TV, < V < — . 

- ' 2 ~ ~ 2 



Remark 2.1 Thus we have harmonized the scales along the fibre lines. 

Definition 2.2 The distance d(P\, P2) between the points P± and P2 is defined by 
the arc length of the shortest geodesic curve from Pi to P2. 

Definition 2.3 The geodesic sphere of radius p (denoted by Sp^p)) with centre 
at the point Pi is defined as the set of all points P2 in the space with the condi- 
tion d(Pi,P 2 ) = P- Moreover, we require that the geodesic sphere is a simply 
connected surface without selfintersection in S 2 xR space. 

Remark 2.4 We shall see that this last condition depends on radius p. 

Definition 2.5 The body of the geodesic sphere of centre Pi and of radius p in 
the S 2 x R space is called geodesic ball, denoted by Bp 1 (p), i.e. Q G Bp 1 (p) iff 
0<d{P u Q)<p. 

In [fT2l we have proved that S(p) is a simply connected surface in E 3 if and only 
if p E [0, 7r), because if p > n then there is at least one v E [— |, |] so that 
y(r,v) = z(t,v) = 0, i.e. selfintersection would occur (see (2.4)). Thus we 
obtain the following 

Proposition 2.6 The geodesic sphere and ball of radius p exists in the S 2 x R 
space if and only if p E [0, tt}. 

We have obtained (see [12]) the volume formula of the geodesic ball B(p) of 
radius p by the metric tensor and by the Jacobian of (2.4): 
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Theorem 2.7 



voi(B( P )) = I — - 1 2 \3/2 dx d y dz = 

fP ri r n 

= J J J \t ■ sin(cos(f )r)| du dv dr = (2.5) 



rp 

2n I I It • sin(cos(u)r)| dv dr. 



o 
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2.1 On fundamental domains 

A type of the fundamental domain of a studied space group can be combined as 
a fundamental domain of the corresponding spherical group with a part of a real 
line segment. This domain is called S 2 x R prism (see |[T2"l0 . This notion will be 
important to compute the volume of the Dirichlet-Voronoi cell of a given space 
group because their volumes are equal and the volume of a S 2 x R prism can be 
calculated by Theorem 2.8. 




Figure 1: Prism-like fundamental domains 



The p-gonal faces of a prism called cover-faces, and the other faces are the 
side-faces. The midpoints of the side edges form a "spherical plane" denoted by 
II. It can be assumed that the plane n is the base plane: in our coordinate system 
(see (2.2)) the fibre coordinate t = 0. 

From [fi2l we recall 

Theorem 2.8 The volume of a S 2 x R trigonal prism Vb b 1 b 2 c c 1 c 2 an d of a dig- 
onal prism Vb BiC Ci in S 2 xR (see Fig. l.a-b) can be computed by the following 
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formula: 

Vol{V) = A-h (2.5) 

where A is the area of the spherical triangle AqAiA 2 or digon A Ai in the base 
plane II with fibre coordinate t = 0, and h = B Cq is the height of the prism. 



3 Ball packings 

By remark (1.2) a S 2 x R space group T has a compact fundamental domain. 
Usually the shape of the fundamental domain of a group of S 2 is not determined 
uniquely but the area of the domain is finite and unique by its combinatorial mea- 
sure. Thus the shape of the fundamental domain of a cry stallo graphic group of 
S 2 x R is not unique as well. 

In the following let T be a fixed glide reflection space group of S 2 x R. We 
will denote by d(X, Y) the distance of two points X, Y by definition (2.2). 

Definition 3.1 We say that the point set 

V(K) = {X G S 2 xR : d(K,X) < d(K s , X) for all g G T} 

is the Dirichlet-Voronoi cell (D-V cell) to T around the kernel point K e S 2 x R. 

Definition 3.2 We say that 

T x = {g G T : X g = X} 
is the stabilizer subgroup of X G S 2 x R in V. 

Definition 3.3 Assume that the stabilizer Yk = I i-e. V acts simply transitively 
on the orbit of a point K. Then let Bk denote the greatest ball of centre K inside 
the D-V cell V[K), moreover let p(K) denote the radius ofBx- It is easy to see 
that 

p{K) = mir^-d(K,K s ). 

The T-images of B K form a ball packing B V K with centre points K G . 
Definition 3.4 The density of ball packing B V K is 

6{K) = VoWiK)- 
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It is clear that the orbit K r and the ball packing B V K have the same symmetry 
group, moreover this group contains the starting cry stallo graphic group T: 

SymK r = SymB r K > T. 



Definition 3.5 We say that the orbit K r and the ball packing B V K is characteristic 
if SymK v = V, else the orbit is not characteristic. 

3.1 Simply transitive ball packings 

Our problem is to find a point K e S 2 xR and the orbit K r for T such that Y K = I 
and the density 5{K) of the corresponding ball packing B V (K) is maximal. In this 
case the ball packing B V (K) is said to be optimal. 

The lattice of V has a free parameter p(T). Then we have to find the densest 
ball packing on K for fixed p(T), and vary p to get the optimal ball packing. 

6(T) = max (8(K)) (3.1) 
k, P (r) 

Let T be a fixed glide reflection group. The stabiliser of K is trivial i.e. we 
are looking the optimal kernel point in a 3-dimensional region, inside of a fun- 
damental domain of T with free fibre parameter p(T). It can be assumed by the 
homogeneity of S 2 x R, that the fibre coordinate of the center of the optimal ball 
is zero. 



3.2 Optimal ball packing to space group 12. 1. 3 

Now we consider the following point group: 

(+,0;[]; {(2,3,4)}) xl R ; 

T := {gl, g2, g 3 - g?, g2, g3> (glg2) 2 , (gigs) 3 , (g2g 3 ) 4 }- 

This is the full isometry group of the usual cube surface, generated by the three 
reflections gj, i = 1,2,3. The possible translation parts of the generators of 
T will be determined by (1.2) and by the defining relations of the point group. 
Finally, from the so-called Frobenius congruence relations we obtain the four non 
equivariant solutions: 

(ri, r 2 , r s ) = (0, 0, 0), (0, 1, 0), (J, 0, 1), (1 1 J). 
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If (ti, r 2, T3) = (|, 0, \) then we get the S 2 x R space group 12. 1. 3. The funda- 
mental domain of its point group is a spherical triangle A AiA 2 with angles | , | , 
I lying in the base plane II (see Fig. 2). It can be assumed by the homogeneity of 
S 2 x R, that the fibre coordinate of the center of the optimal ball is zero and it is 
an interior point of A AiA 2 triangle. 

We shall apply the Cartesian homogeneous coordinate system introduced in 
Section 2 (see Fig. 2) and the usual geographical coordiantes (0, 9), (— n < 4> < 
7r, — ~ < 9 < ~ ) of the sphere with the fibre coordinate t G R (see (2.2)). 

We consider an arbitrary interior point K{ ) of spherical triangle 

A AiA 2 in the above coordinate system in our model by the following equations: 

x° = 1, x 1 = cos cos 6* , x 2 = sin0cos#, x 3 = sin# (3.3) 

Let B r (R) denote a geodesic ball packing of S 2 x R space with balls B(R) of 

♦ z 




Figure 2: 



radius R where their centres give rise to the orbit K . In the following we consider 
to each ball packing the possible smallest translation part t(K, R) (see Fig. 2) 
depending on T, K and R. A fundamental domain of Y is its D-V cell T>(K) 
around the kernel point K. It is clear that the optimal ball Bk has to touch some 
faces of its D-V cell. The volume of D(K) is equal to the volume of the prism 
which is given by the fundamental domain of the point group T of T and by the 
height 2t(R, K). The images of D(K) by our discrete isometry group covers the 
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-o.£ 



Figure 3: 

S 2 xR space without overlap. For the density of the packing it is sufficient to relate 
the volume of the optimal ball to that of the solid D(K) (see Definition 3.4). 

It is clear, that the densest ball arrangement Br(R) of balls B(R) has to hold 
the following requirements: 

(a) d(K, K 92 ) = 2R = d(K, K T91 ), 

(b) d(K, K 92 ) = 2R = d(K, K T93 ) } 

(c) d(K, K 2r ) > 2R (3.5) 
(d) Balls of radius R with centres 
K, K 92 , K T9 \ K T9 \ K 2t form a packing. 

Here d is the distance function in the S 2 x R space (see Definition 2.2). The 
equations (a) and (b) mean that the ball centres K T91 and K T93 lie on the equidis- 
tant geodesic surface of the points K and K 2r which is a sphere in our model in 
this case (see 0). 

We consider two main ball arrengements: 

1. We denote by Br(Ro, K ) those packing where requirements (3.5) and 
d(K, K 2t ) = 2R hold (see Fig. 3). 

2. We denote by £> r (i? 1; Ki) those packing where requirements (3.5) and 
d(K T9 \K T93 ) = 2R hold (see Fig. 4). 
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Figure 4: 

First we determine the coordinates of the points Ki, (i = 1, 2) (K, t is given by 
(3.3) with parameters 4> and 9), the radius R of the ball, the volume of a ball B(R) 
and the density of the packing in both main cases. We get the following solutions 
by systematic approximation, where the computations were carried out by Maple 
V Release 10 up to 30 decimals: 

00 ~ 0.24389626, 9 « 0.20663860, R « 0.23860571, 
VoI(B(Rq) « 0.05668684, 5(R ,K ) « 0.45373556. LO) 

01 w 0.30773985, 6 l w 0.17313169, i2 x w 0.30299179, 
Vol{B(R 1 ) RJ 0.11580359, 5(^1,^) w 0.44472930. ( "'' 6) 

We obtain by careful investigation of the density function i^) (i? e [i?o, -Ri]) 
of the considered ball packing the following: 

Theorem 3.6 The ball arrangement Br(Ro, K ) (see Fig. 3) provides the densest 
symply transitive ball packing belonging to the S 2 x R space group 12. 1. 3. 

3.3 The densest simply transitive ball packing 

We consider the following point group: 

(+, 0, [ ]{(?,?)})xl R , q>2; 

To = (gi,g 2 -gi,gi,(gig 2 ) 9 )- 



12 



Jend Szirmai 




Figure 5: 



This point group is generated by two reflections gj, i = 1,2,3. The possible 
translation parts of the generators of T will be determined by (1.2) and by the 
defining relations of the point group. Finally, we obtain from the so-called Frobe- 
nius congruence relations three non equivariant solutions: 

(r 1; r 2 ) = (0,0), (0,i), {\\). 

If (ti, r 2 ) = (|, |) then we have obtained the S 2 x R space group 2q. I. 2. 

The fundamental domain of the point group of the considered space group is 
a spherical digon A Ai with angle - in the base plane IT. Similarly to the above 
section can be assumed, that the fibre coordinate of the center of the optimal ball 
is zero and it is an interior point of A Ai digon (see Fig. 5). 

In the following we consider ball packings belonging to q = 2. We use also 
the above introduced Cartesian homogeneous coordinate system and the usual 
geographical coordinates (</>, 9), (— n < <p < it, — | < 9 < | ) of the sphere with 
the fibre coordinate t G R (see (2.2)). 

We consider an arbitrary interior point K(l, x 1 , x 2 , x 3 ) = K((p, 9) of spherical 
digon AqAi in the above coordinate system in our model (see Fig. 5). 

Our aim is to determine the maximal radius R of the balls, and the maximal 
density 5(K, R). 

The ball arrangement B opt (K, R) is defined by the following equations: 



(a) d(K, K T91 ) = 2R = d(K, K T92 ), 
(6) d{K T9l ,K T92 ) = 2R, 



(3.7) 
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We can determine the coordinates of the point K, the radius R of the ball, the 
volume of a ball B(R) and the density of this packing: 




Figure 6: The densest simply transitive geodesic ball packing 



eh =- « 0.78539816, 6 = 0, R =-« 1.57079633, 

V 4 2 ' (3.8) 

Vol{B(R) « 13.74539472, if) « 0.80407553. 

Similarly to the above Section we can prove the following theorem: 

Theorem 3.7 The ball arrangement B opt (R, K) (see Fig. 6) provides the densest 
simply transitive ball packing belonging to the S 2 x R space group 2q. I. 2. 

Considering all in this paper investigated space groups and in lfi~2ll discussed gen- 
eralized Coxeter S 2 x R space groups we get the next theorem: 

Theorem 3.8 The ball arrangement B opt (R, K) provides the densest simply tran- 
sitive ball packing belonging to the generalized Coxeter and glide reflections gen- 
erated S 2 x R space groups. 

By Theorems 2.7 and 2.8 and by Definitions 3.3 and 3.4, similarly to the above 
space groups, we have determined the data (radii, densities and volumes of opti- 
mal balls) of the optimal simply transitive ball packings to each glide reflections 
generated S 2 x R space group which are summarized in Table 1. 
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Remark 3.9 The space groups 2q. I. 2, 2qe. I. 3, 4q. I. 2, 4q. I. 3, 4q. I. 4, 

4qe. I. 5, 4qe. I. 6 depend on parameter q thus their optimal ball packings depend 
also on q but in the Table 1 we give only the data of the densest ball packing 
indicating its q parameter to each considered space group. 



Table 1 



Space group 


R Vol(B K (R)) 5 


2q. I. 2, q — 2 


f ps 1.57079633 ps 13.74539472 ps 0.80407553 


2qe. I. 3, q — 2 


f ps 1.57079633 


w 13.74539472 


ps 0.69634983 


4q. I. 2, q — 2 


ps 0.64360446 


w 1.08624788 


ps 0.53722971 


4q. I. 3, q — 2 


ps 0.67517586 


w 1.25058159 


ps 0.58958340 


4q. I. 4, q — 2 


ps 0.95531662 


w 3.43551438 


ps 0.74837055 


4qe. I. 5, q = 2 


ps 0.64360446 


w 1.08624788 


ps 0.53722971 


4qe. I. 6, q = 2 


ps 0.67517586 


w 1.25058159 


ps 0.58958340 


11. I. 2 


ps 0.46364761 


« 0.41154972 


ps 0.58861600 


12. I. 2 


« 0.22770028 


w 0.04928081 


ps 0.41334779 


12. I. 3 


ps 0.23860571 


« 0.05668684 


ps 0.45373556 


12. I. 4 


w 0.31004511 


« 0.12404486 


ps 0.53597559 


13. I. 2 


« 0.18705243 


ps 0.02735051 


ps 0.49222087 



It is timely to arising the above question for further space groups in S 2 x R 
space. 

In this paper we have mentioned only some problems in discrete geometry of 
S 2 xR space, but we hope that from these it can be seen that our projective method 
suits to study and solve similar problems ([j7), IfTTl . Ifi~2l ). Analogous questions in 
other homogeneous Thurston geometries are interesting. 
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